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We have been studying effects of dynamical quarks on various hadronic observables, using our recently formu- 
lated improvement for staggered fermions. To illustrate improvement, we show that the light hadron spectrum 
in the quenched approximation gives remarkably good scaling. We highlight three new results: (1) We find no 
apparent quark loop effects in the Edinburgh plot with 2+1 flavors of dynamical quarks at a — 0.14 fm. (2) We 
show that dynamical quarks modify the shape of the heavy quark potential. (3) We present results hinting at 
meson decay effects in light hadron spectroscopy. 



1. INTRODUCTION 

Recent proposals for Symanzik improvement of 
the staggered fermion action have produced en- 
couraging results, particularly in an improved fla- 
vor and rotational symmetry (l]j|] . Adding a few 
terms to the conventional Kogut-Susskind action, 
namely three-link, five- link and seven-link staples 
and a third-neighbor coupling, removes all tree- 
level 0(a 2 ) errors The computational price 
for using this action is quark mass dependent, and 
is roughly a factor of 2.5 for the lightest masses 
we are using. Here we extend previous work with 
our preferred "Asqtad" action, studying scaling in 
the quenched approximation, and exploring the 
effects of dynamical quarks on the light hadron 
spectrum and the heavy-quark potential. 

We have accumulated a library of gauge config- 
urations, both quenched (one-loop Symanzik im- 
proved gauge action) and dynamical.^] Our 
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1 The MILC code, including the improved KS action code, 
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dynamical fermion lattices are generated in the 
presence of two lighter flavors (u, d) and one heav- 
ier flavor (s). For the dynamical fermion sample 
the gauge coupling is tuned so that the lattice 

is publicly available. Gauge configurations are also avail- 
able. Please contact doug@physics.arizona.edu or de- 
tar@physics.utah.edu for details 
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Figure 1. Flavor symmetry test. Mass of 
the non-Goldstone local pion vs lattice spacing 
in units of r\ (see text) for progressively im- 
proved actions: (octagon) conventional staggered 
fermion and single-plaquette gauge action, (dia- 
mond) same but with an improved gauge action, 
(cross) same but with the Naik third-neighbor 
term, and (square) our preferred Asqtad action. 

spacing remains the same as the quark mass is 
varied, allowing an exploration of quark mass ef- 
fects independent of scale. The current parameter 
set is shown in the table above. 

2. SCALING TESTS 

To measure the degree of residual flavor sym- 
metry breaking, we compute the mass of the local 
non-Goldstone pion (7105) as a function of lattice 
spacing at a fixed value of the Goldstone pion 
mass {vtlq) on a series of quenched lattices |TT| . 
In particular, we fix mari = 0.807(3) (Goldstone 
7r) with the scale r\ set by the force -pQQstatic be- 
tween static quarks T , iF , qQ static (ri) = 1, a variant 
of Sommer's scale[[[o| . The value of ro / r% depends 
on the quark masses, varying from 1.376(2) for 
the quenched theory to 1.44(1) for physical quark 
masses. The effect of improving the fermion ac- 
tion is shown in Fig. [|. It is expected that the 
curves for all actions extrapolate to the Goldstone 



Figure 2. Scaling test. Nucleon mass vs lattice 
spacing. Same as [l], but including Wilson-clover 
results. 
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Figure 3. Same as|^, but for the rho meson mass, 
pion mass. 

Scaling of the p and nucleon masses is shown in 
Figs. H and |[ Results for variants of the Wilson- 
clover action are also shown jl2|,[l3]]. For these 
masses, the Asqtad action gives the best result. 
Nonetheless, it is encouraging that staggered 
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Figure 4. Edinburgh plot comparing results for a 
quenched and 2 + 1 flavor dynamical simulation 
at a = 0.14 fm. Octagons show the experimental 
value (lower left) and infinite quark mass value 
(upper right). 

fermion and Wilson-clover simulations agree in 
the continuum limit. 

3. DYNAMICAL QUARK EFFECTS: 
EDINBURGH PLOT 

To determine the effect of dynamical quarks 
on the light hadron spectrum, we ran a series of 
quenched and 2+1 dynamical fermion simulations 
with gauge couplings tuned to fix the lattice spac- 
ing a — 0.14 fm, using 20 3 x 64 lattices for both 
the quenched and dynamical runs. The plotted 
points arc from the runs described in Table I with 
light quark masses down to 0.02a _1 and from the 
quenched run at 10/.g 2 = 8.0. The resulting Ed- 
inburgh plot in Fig. ^ shows no discernible change 
when dynamical quarks are introduced. This 2+1 
flavor a = 0.14 fm result appears to be in conflict 
with our previous claims, based on the conven- 
tional action with two sea quark flavors in the 
a — > limit ||l6| . We are currently investigating 
whether the difference comes from the number 
of flavors, the continuum extrapolation, or some 
other source of systematic error. 



Figure 5. Square root of the string tension in 
units of r\ vs the squared ratio of the pion to rho 
mass. The octagons are done with three degener- 
ate flavors (for m q > m s ), while the squares are 
runs with two light flavors and a fixed strange 
quark mass. 

4. DYNAMICAL QUARK EFFECTS: 
HEAVY QUARK POTENTIAL 

Dynamical quark loops modify the heavy quark 
potential flFfl] , e.g., by opening the two-meson de- 
cay channel and by altering the running of the 
Coulomb coupling. In Fig. |^ we plot the heavy 
quark potential with and without 3 flavors of sea 
quarks of mass am q ~ 0.05. The vertical scale is 
adjusted to give agreement at r±, and since the 
force at r\ is used to set the scale, the two poten- 
tials have the same slope there. 

To the extent the potential changes shape un- 
der the influence of sea quarks, setting the lattice 
scale from the heavy quark potential clearly has 
limitations. This point is illustrated in Fig. |^, 
where we plot the ratio of two possible heavy- 
quark-potential-inspired scales as a function of 
the ratio of the pion to rho mass. 

5. DYNAMICAL QUARK EFFECTS: 
MESON DECAY 

As the pion mass drops to zero with decreas- 
ing quark mass, one expects to see evidence for 
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Figure 6. Heavy quark potential with and without 2 
physical units for the two matched simulations. 

the decay of heavier mesons, such as p — > 2tt 
and a — ► 7777. The mixing between a resonant 
state and a two-meson decay state should result in 
avoided level crossing in the spectrum of the res- 
onant state. While such decays are natural when 
dynamical quarks are present, quenched simula- 
tions may also show hints of decay in some cases, 
through "hairpin" diagrams, although at different 
rates. 

Previous attempts to observe effects of meson 
decay in the p spectrum were unsuccessful fl8| . 
Here we report results of a preliminary attempt 
to observe effects of the decay ao — > rjir. This 
threshold is easier to approach than the p — > tt+tt 
threshold because the final pseudoscalars don't 
have to carry away angular momentum, so they 



1 flavors of sea quarks. Rulers show lattice and 



can both be at zero momentum. However, it is 
difficult, because with staggered quarks the ao ap- 
pears in the same propagators as the pions, so we 
must extract the alternating exponential from a 
propagator containing a much larger simple ex- 
ponential. We measure the ao mass (i.e. lowest 
energy in the ao (0 ++ ) channel) as a function of 
quark mass. This is done on quenched and dy- 
namical fcrmion configurations. Preliminary re- 
sults, shown in Fig. [?], suggest that the lowest ++ 
mass drops more rapidly with decreasing quark 
mass in the dynamical fermion simulation. The 
drop appears to follow the threshold energy of the 
decay channel, also shown. Such a result would 
be expected from avoided level crossing. In fact, 
for the lighter quark masses it is very difficult to 
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Figure 7. Mass of the ++ state seen in our spec- 
trum analysis vs quark mass with and without 
dynamical quarks. The r] — n threshold is plot- 
ted. 

extract masses for the quenched ag. This may 
be a signal of the non-unitarized (because of the 
quenching) couplings to two particle states. 



6. SUMMARY 

We continue studies using a recently proposed 
improved staggered fermion action. We find that 
the light hadron spectrum in our currently pre- 
ferred Asqtad action scales very well. We have be- 
gun a systematic study of quark loop effects with 
2 + 1 flavors. Our present sample of gauge config- 
urations has been constructed so that the lattice 
spacing is kept constant as the quark masses are 
varied. We see no quark loop effects in the light 
hadron spectrum with 2 + 1 flavors at a lattice 
spacing of a = 0.14 fm. However, sea quarks evi- 
dently modify the heavy quark potential both at 
short and long distance, suggesting caution in us- 
ing this potential to set a precision lattice scale at 
unphysical quark masses. Preliminary determina- 
tions of the mass of the do as a function of quark 
mass suggest evidence for the decay to rj + ir. 
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